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We have carried out extensive first principles doping-dependent computations of angle- resolved 
photoemission (ARPES) intensities in L^-xSr^CuCU (LSCO) over a wide range of binding ener- 
gies. Intercell hopping and the associated 3-dimensionality, which is usually neglected in discussing 
cuprate physics, is shown to play a key role in shaping the ARPES spectra. Despite the obvious 
importance of strong coupling effects (e.g. the presence of a lower Hubbard band coexisting with 
mid-gap states in the doped insulator), we show that a number of salient features of the experi- 
mental ARPES spectra are captured to a surprisingly large extent when effects of fc z - dispersion are 
properly included in the analysis. 

PACS numbers: 79.60.-i, 71.18.+y, 74.72.Dn 



La2- x Sr x Cu04 (LSCO) has drawn intense interest as 
a model system for understanding one of the most hotly 
debated issues in condensed matter physics currently, 
namely, how does a Mott insulator, La2CuC>4 (LCO), 
develop into a superconductor when doped with holes 
via La/Sr substitution, and what is the route taken by 
the electronic states in the insulator to achieve this re- 
markable transformation into a metal. 1,2 Angle-resolved 
photoemission spectroscopy (ARPES) has been brought 
to bear on these questions in recent years as techniques 
for preparing LSCO surfaces have improvedi^^iL^ 
ARPES spectra in LCO find a lower Hubbard band 
(LHB) associated with an insulating state, which persists 
with finite hole doping, losing intensity without shifting 
in energy. At the same time new states - the so-called 
mid-gap states - appear close to the Fermi level, which 
evolve into the conventional metallic bands near opti- 
mal doping. Doping-dependent ARPES spectra in LSCO 
thus address a wide range of issues concerning stripe 
and pseudogap physics and their relation to the mech- 
anism of high-temperature superconductivity, and bear 
on questions of non-Fermi liquid behavior or gossamer 
superconductivity^*^, among others^i 

We have recently shown with the example of Bi2212 
that the residual /^-dispersion of bands in a quasi- 
2D material will induce an irreducible linewidth in 
ARPES peaks, which is unrelated to any scattering 
mechanisms J" 2 " This effect becomes accentuated in LSCO, 
where the bands possess a greater 3D character compared 
to Bi2212. This article reports extensive first principles 
computations of the ARPES intensity in LSCO with the 
goal of ascertaining the extent to which /^-dispersion af- 
fects the experimental spectra of the mid-gap and LHB 
states. The calculations properly model the photoemis- 
sion process and include the crystal wave functions to 
describe the initial and final states in the presence of 
the surface and take account of the associated dipole 
matrix element and its dependencies on photon energy 
and polarization. The interplay between the effects of 
the ARPES matrix elemen1si2*i2ii£ and the fc z -dispersion 
yields computed intensity maps for emission from the 



Fermi energy (Ep) as well as for binding energies sev- 
eral hundred meV's below Ep, which are in surprising 
accord with the corresponding measurements. 

Our calculations give insight into a number of salient 
features of the experimental spectra such as, the disper- 
sion of the mid-gap states and the characteristic broad- 
cnings and symmetries or lack thereof in the photointen- 
sities. Evidence of physics beyond the framework of the 
conventional local density approximation (LDA) based 
picture is especially clear in the strongly underdoped 
regime, both in the appearence of a d-wave-like pseu- 
dogap in the mid-gap states and in the presence of the 
LHB near half filling. We discuss this insulating regime 
in terms of tight-binding (TB) computations, in which 
the LDA-inspired overlap parameters are supplemented 
with a Hubbard U or a pseudogap A. The effect of fcz- 
dispersion is included through an intercell hopping pa- 
rameter - to our knowledge for the first time in connec- 
tion with a TB description of the cuprates. In this way, 
via comparisons between theory and experiment, we ad- 
duce a strong connection between the LDA generated 
metallic states and the evolution of the mid-gap band 
throughout the doping range. 

With regard to relevant computational details, the 
fully self-consistent electronic structure of tetragonal 
LCO was obtained within the LDA by using the well- 
established all electron Green function methodology 16 ; 
our band structure and Fermi surface (FS) are in accord 
with published data 1 - . Following common practice, the 
metallic state of LSCO obtained by doping LCO with Sr 
is assumed to be described by the LDA generated metal- 
lic band structure for LCO. The effect of La/Sr substitu- 
tion is mainly to adjust the electron concentration and, 
therefore, we have modeled LSCO at any given cc-value 
by invoking a rigid band filling of the band structure of 
LCO with the appropriate number of electrons per unit 
cell. All presented ARPES intensities have been com- 
puted within the one-step photoemission formalis m, as- 
suming an LaO-layer terminated surface; see Refs. Ibll ll 
and Il2l for details. 

Fig. 1 considers emission from Ep for optimally doped 
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FIG. 1: (Color) ARPES spectra for emission from the Fermi 
energy in optimally doped LSCO for x — 0.15. (a): Red filled 
region on left hand side gives the 3D Fermi surface projected 
onto the (k x , fc a )-plane, and denotes the area of allowed emis- 
sions. The right hand side gives computed intensity including 
the effect of the ARPES matrix elements with red lines mark- 
ing boundaries of the red filled regions, (b): Corresponding 
experimental spectral. 



LSCO and shows the remarkable degree to which the FS 
map can be understood within the framework of the con- 
ventional LDA-based picture. We start by looking at the 
left side of (a), where the filled red region gives the pro- 
jection of the FS onto the (k x , fc ?/ )-plane and encompasses 
various FS cross-sections as a function of k z ^ Emission 
of photoelectrons is possible in principle from any part of 
this red region due to the quasi-2D nature of states. [In 
a strictly 2D system, the red region will collapse into a 
standard FS contour of zero width.] We emphasize that 
the photointensity within this red region will in general 
not be uniform as it will be modified by the effect of 
the ARPES matrix element. This aspect is delineated 
on the right side of (a), where the theoretically predicted 
ARPES intensity is shown superposed with the bound- 
aries of the region of allowed transitions by red lines. 

When we compare the theoretical intensities in 
Fig. 1(a) with the corresponding experimental results^i^ 
in Fig. 1(b), the most striking feature is the appearance of 
wing-like structures around the anti-nodal points M(tt, 0) 
in both theory and experiment. On the other hand, along 
the nodal direction in the first Brillouin zone (BZ), T to 
X(t:, 7t), the high computed intensity is not reproduced 
in the measurements. However, we find that this com- 
puted nodal intensity is quite sensitive to photon energy. 
It will be necessary to carry out ARPES experiments over 
a range of photon energies in order to ascertain whether 
or not absence of nodal intensity represents a significant 
effect of correlations beyond the LDA in LSCO. 

Figure 2 expands the preceding discussion to include 
the underdoped and overdoped regimes. The computa- 
tions in all cases simulate experimental conditions^ of 
resolution, polarization of light and photon energy. The 
theory is seen to provide a good overall description of the 
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FIG. 2: (Color) Similar to Fig. 1, except that this figure 
compares theoretical and experimental ARPES intensities in 
the underdoped x = 0.06 [panels (a) and (b) (Ref. 0)] and 
overdoped x — 0.22 [panels (c) and (d) (Ref. 5j)] cases. 

data over the entire doping range, discrepancies along 
the nodal direction notwithstanding. Around the antin- 
odal point, incorporation of fc z -dispersion allows the evo- 
lution of the FS across the (broadened) Van Hove sin- 
gularity (VHS) to be analyzed in detail. In the under- 
doped sample the intensity is concentrated in two fea- 
tures lying above and below the M(tt, 0)-point. In con- 
trast, the FS crosses the VHS near optimal doping, and 
the overdoped sample displays a greater spectral weight 
around the M(ir, 0)-point. Since the midgap states evolve 
into conventional bands with increased doping, the good 
agreement with theory is perhaps to be expected in op- 
timally and overdoped samples, but the continued agree- 
ment for the underdoped sample is quite remarkable - 
particularly since the midgap states have a pseudogap, 
as discussed below. 

Figs. 1 and 2 make it clear that the broadening of 
ARPES spectra resulting from the effect of /c z -dispersion 
is essentially zero along the nodal direction and that it 
increases only gradually as one moves towards the antin- 
odal region. By contrast in the antinodal region, the ef- 
fects of the VHS conspire with those of the /^-dispersion 
to produce effective linewidths which increase rapidly as 
one moves away from the antinodal point. The antin- 
odal point itself is anomalous in that the broadening can 
be quite small parallel to the BZ boundary (X(ir, it) to 
M(tt, 0)) as seen for example from Fig. 2(a). The present 
analysis shows that, even if correlation effects were ab- 
sent, features in the antinodal region would generally be 
considerably broader than nodal ones, due simply to the 
effect of fc z -dispersion. 

Caution should be exercised in interpreting the sig- 
nificance of differences between theory and experiment 
in Figs. 1 and 2 around say the M(ir, 0)-point or for 
that matter even the nodal direction. On the theoret- 
ical side, some of these details are sensitive to photon 
energy and to the precise position of the Fermi energy in 
relation to the VHS in the density of states. Similarly, 
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FIG. 3: (Color) Band dispersion in optimally doped LSCO, 
x — 0.15. (a): Computed ARPES intensities along the F — 
M — X — r -line in the BZ; color scale is the same as that 
in Fig. 1. Red lines give bands corresponding to different 
values of k z : k z = (solid line), k z = tt/c (dotted), and 
k z = 2ir/c (dashed), (b): Second derivative of the theoretical 
spectra in (a) as a grey scale plot, (c): Experimental second 
derivatives- corresponding to the theoretical results in (b). 
Note scale change between panels (b) and (c). 
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FIG. 4: (Color) Theoretical band structures of LSCO are 
shown overlayed on experimental ARPES spectra (second 
derivatives) of Ref. for two different dopings: (a) x = 0.03, 
and (b) x = 0.0. The data are given on a hot color scale in 
which red and black denote lows. Two sets of tight-binding 
bands are shown by green and white lines as detailed in the 
text. The two different bands in each case correspond to two 
values of k z : k z — and k z = 2n/c. Different line types de- 
note computed relative spectral weights varying from 1.0-0.8 
(solid lines), 0.8-0.2 (dashed) and 0.2-0.01 (dotted). 



higher resolution experiments are needed to pin down 
some of the fine details in the experimental spectral 
Ultimately, ARPES evidence for exotic physical effects 
beyond the conventional LDA-type pictures (e.g. stripes, 
pseudogaps, strong correlations) must derive from a care- 
ful analysis of the differences in the fine structure between 
theoretical and experimental intensities such as those of 
Figs. 1 and 2. 

The typical behavior of spectra below Ep is consid- 
ered in Fig. 3 with the example of the optimally doped 
system. Much of the commentary in connection with 
Figs. 1 and 2 above is applicable and need not be re- 
peated. In Fig. 3(a), the computed ARPES intensity 
along the T — M as well as the T — A"-directions lies 
within the limits given by the bands for various fc z -values, 
albeit with a strong modulation by the ARPES matrix 
element. Once again, the spectral lines are substantially 
broader along T — M compared to T — X. Theoretical 
and experimental^ results in (b) and (c), respectively, are 
in substantial accord in this respect. Note differences in 
vertical scale in (b) and (c) , indicating a renormalization 
(i.e. reduction) in bandwidth by roughly a factor of two 
over the LDA predictions, which is quite common in the 
cuprates. Also, the measured ARPES intensity in (c) is 
seen to pull away from the Fermi energy, which is not 
the case in the theoretical plot of (b). This difference 
however is to be expected since the theory refers to the 
normal state, while the measurements are taken from a 
superconducting sample. 

We have fitted the first principles FS and band disper- 
sions near Ep in LSCO within the TB model using the 



formic 

e k = —2t(c x + Cy) - At 'c x Cy — 2t'' '(c 2x + c 2y ) 

-2 T z (k\\) c z (c x - c y ) 2 , (1) 

with Ci = cos(kia) and en — cos(2kia), i = x,y, c z = 
cos(k z c/2), and 

T z = t z cos{k x a/2) cos(k y a/2). (2) 

The extra angular dependence in Eq. 2 accounts for the 
staggered stacking of neighboring Cu02-planes, and re- 
sults in vanishing fc z -dispersion along the BZ boundaries. 
The parameters which yield a good fit are: t = 0.32eV, 
t' = -0.25t, t" = -0.02i, and t z = 0.16i. The lat- 
ter value should be compared with t z = O.lli estimated 
from transport In order to simulate the presence of 
a pseudogap and/or the superconducting gap, the bare 
bands of Eq. 1 should be replaced by: e k — > ±E k , where 
E% — e l + and A k — A Q (c x - c y )/2. In this way 
a d-wave-like gap with a maximum value of Ao at the 
M(tt, 0)-point is produced in the spectrum; we find that 
Ao grows with underdoping as a pseudogap. 

Fig. 4 considers the underdoped regime below the 
superconductor-insulator (SI) transition (around x — 
0.06). Focusing on the experimental^ spectra first, the 
half-filled case (x = 0.0) in (b) is seen to display the 
presence of the lower Hubbard band (LHB) at a binding 
energy of « 0.5 eV, and little weight around Ep. When 
the system is doped with a small amount of holes, we 
see from (a) the appearance of a substantial weight in 
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the mid-gap spectrum near the M(%, 0) and (ir/2,ir/2) 
points, while the LHB is essentially unchanged. 

The experimental results of Fig. 4 cannot of course be 
understood within the conventional LDA-based picture, 
which fails to produce the insulating state at half-filling. 
Insight can however be obtained via the TB model, and 
accordingly, we have carried out two sets of such calcula- 
tions: (i) The LHB is fit to TB mean-field computations 
for a saturated antiferromagnet with Hubbard parameter 
U = Qt, with all other parameters taken directly from the 
LDA bands. The fcz-dispersion is accounted for via the 
t z parameter (see Eqs. 1 and 2), which to our knowledge 
has not been included in any previously published work. 
These TB bands for the two extremal fc z -values are given 
by green lines in Fig. 4, with different line types giv- 
ing associated spectral weights (see figure caption) . Our 
parameter values are consistent with earlier data from in- 
sulating cuprates based on the ARPESS and spin wave 
spectra 2 ^, (ii) TB calculations with the same parame- 
ters as those used in the preceding case, except that the 
Hubbard parameter U is replaced by a d-wave-like pseu- 
dogap with Ao = 220 meV at M(tt, 0), which accounts 
for the presence of such a gap in the experimental spectra 
at x = 0.03. 

The theoretical bands of Fig. 4 provide a handle on 
understanding some salient features of the experimental 
spectra. We consider the LHB along the T — M and 
M — X lines first. At both x — and x = 0.03, the 
measured LHB is seen to extend over the binding energy 
range of 0.4-0.8 eV, more or less symmetrically about the 
M(tt, 0)-point, at least insofar as its width is concerned. 
Much of this width can be ascribed to the effect of k z - 
dispersion, which will induce spectral lines to spread be- 
tween the boundaries given by the pair of green bands. 
Note also that the computed LHB is symmetric around 
M(n, 0) due to the effect of zone-folding in the AFM insu- 
lating state. The calculated spectral weights of the green 
bands are seen to be larger along T — M compared to the 
M — X line^ The experimental broadening at M(tt, 0) 
is anomalously large and not due to fc z -dispersion. This 
however is not surprising since the M(jr, 0)-point lies on 
the AFM zone boundary and is susceptible to 'hot spot' 
scattering. The situation along the T — X line is sharply 
different. Here the effect of fc z -dispersion is negligible. 
Moreover, r — X is orthogonal to the AFM zone bound- 
ary. Therefore, substantial observed widths of the LHB 
along r — X cannot be due to the effect of either k z - 
dispersion or of magnetic scattering. 

We turn now to comment on the mid-gap spectrum 
in Fig. 4(a). Here the aforementioned symmetry of the 
LHB with respect to the M(tt, 0)-point is absent. The in- 
tensity of mid-gap states cuts off quite abruptly near the 
M-point along the M — X line. This is in accord with the 
LDA-based bands (white lines), which show little effect 
of fc z -dispersion and also possess little weight to the right 
hand side of M(n, 0). On the other hand, a broad patch 
of intensity is seen around M(tt, 0) extending toward T. 



This is to be expected due to the effect of fc 2 -dispersion 
in view of the boundaries given by the pair of white LDA- 
based bandsi^i Interestingly, along the r — X -line, the 
experimental mid-gap band is quite sharply defined, con- 
sistent with the fact that the LDA bands display little k z 
induced broadening. 

Taken together, the comparisons of Figs. 1-4 paint a 
remarkable picture of the way mid-gap states evolve in 
LSCO with doping. The fact that the spectrum of the 
lightly doped insulator along the T — X line in Fig. 4(a) 
essentially follows the LDA bands, suggests that even the 
first mid-gap states created when holes are added into 
the Mott insulator mimic metallic states^ However, the 
precise nature of the mid-gap states and how the asso- 
ciated pseudogap evolves to yield a Luttinger-like Fermi 
surface remains an important theoretical question. The 
fact that hopping appears only weakly renormalized with 
doping, but that the spectral weight or the intensity of 
the mid-gap band scales with x, is much in the spirit of 
a gossamer type mode&iS. The presence of a pseudogap 
near (jr, 0) will result in predominantly incoherent c-axis 
transport in the underdoped cuprates2&, whereas there 
may be coherent c-axis transport in overdoped samples. 27 
Certainly, the ARPES spectra in the vicinity of Ep are 
described in considerable detail by the LDA computa- 
tions as seen from Figs. 1 and 2. 

In conclusion, we have shown with the example of 
LSCO, that effects of fc z -dispersion, which have been ne- 
glected in most of the existing literature on the cuprates, 
play a key role in explaining the observed broadening 
in the ARPES spectra throughout the BZ for both the 
mid-gap band as well as the LHB. These results provide 
a new benchmark for testing strong correlation models 
of the ARPES spectra. For example, evidence for stripe 
or marginal Fermi liquid physics^ must take into ac- 
count the w-dependent broadening associated with fcz- 
dispersion. Finally, despite the obvious importance of 
strong correlation effects, the remarkable extent to which 
simple LDA-type metallic states describe the dispersion 
of the mid-gap spectrum in LSCO is an observation call- 
ing for theoretical interpretation. Our results will pro- 
vide a deeper understanding of strong correlation effects 
in the cuprates, including insight into the mid-gap states 
and the pseudogap, and will likely require some revisions 
of associated models. 
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